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I. INTRODUCTION 

Modern applications of quantum theory renewed interest in characterization of the set of mixed quantum states. It 
is often necessary to quantify, how a certain mixed state may be approximated by another one. For this purpose one 
may use various distances defined in the space of mixed quantum states, e.g. the trace distance, the Hilbert-Schmidt 
(HS) distance and the Bures distance Ij. The latter_ distance, is a function of fidelity and may be considered as 
a generalization of the overlap between pure states !3| . Fidelity between a pair of two mixed states is equal to unity 
if and only if they do coincide. Optimal fidelity between the input and the output states are used to characterize the 
quality of a cloning algorithm 0|, . 

In various quantum problems one computes the mean fidelity, averaging over a certain ensemble of states. Analyzing 
dynamics of quantum chaotic systems one studies the fidelity between a given initial state and its image under the 
action of the quantum map ^ . To extract global information on the entire system one averages this quantity over the 
natural, Fubini-Study (FS) measure on the space of initial pure states. Mean fidelity between two images of a random 
pure state obtained by the action of a certain unitary operation and a certain irreversible, contracting superoperator 
£ was found by Bowdery et al. for qubits and by Nielsen 3 for the more general case of TV dimensional states. 
Another expression for the average gate fidelity was recently derived by Bagan et al. jQi] , in which the nontrivial term 
was shown to be proportional to the mean value of the Hilbert-Schmidt scalar product (Tr7i£(7i))i of the SU{N) 
generators 7^ and their images averaged over all of them, i = 1, . . . , N'^ — 1. 

In this work we study a complementary problem and analyze mean fidelity {F) between two independent, random 
states. Such a problem may be interesting from a practical point of view: The average fidelity may serve as reference 
value by analyzing the mean (maximal/minimal) fidelity achieved in a certain protocol of quantum cloning. For 
instance, if the measured (computed) fidelity between a certain pair of one-qubit mixed states only slightly exceeds 
the average value, one should not conclude that these states are more correlated than two generic, random mixed 
states. 

In this way we are in position to propose a general tool measuring the quality of a given theoretical scheme of 
quantum information processing or its experimental realization. Let F denote the mean fidelity between the state 
obtained in an analyzed cloning protocol and the target state. Then the quality of the protocol may be gauged by a 
dimensionless coefficient a — [F — (F)]/ ^ (F^) — (F)^, where {F) denotes the average over an appropriate ensemble 
of random density matrices. 

We are also interested in the probability distribution P{F). Two averaging schemes should be distinguished. In the 
symmetric case, both states are drawn at random according to the same measure (e.g. both pure states or both mixed 
states generated with respect to the same probability distribution). In the non-symmetric case, both probability 
distributions are different: in particular we study the mean fidelity between a random pure state and a random mixed 
state, generated with respect to a specified measure. 

Let us emphasize in this place that there exists no single, naturally distinguished probability measure in the set of 
mixed quantum states [l^, UM, ■ Guessing a mixed state of size N at random we might use additional information, 
if available. For instance, knowing that the mixed state has arisen by the partial tracing over a K dimensional 
environment, the induced measure 11] should be used. In particular, if the size of the environment K and the size 
N of the system are equal, then the measure induced by partial trace coincides with the Hilbert-Schmidt measure. 
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Without any information concerning the random state whatsoever, it will be legitimate to make use of the Bures 
measure, related to the Jeffrey's prior, statistical distance and the distinguishability. 

The paper is organized as follows. In section 2 we review the basic properties of fidelity, while in section 3 we 
introduce the necessary measures in the set of mixed quantum states of size N. Mean fidelity between pure and 
mixed states are analyzed in section 4. Main results of this work are presented in section 5, in which we compute the 
mean fidelity, averaged over two independent random mixed states, generated with respect to an arbitrary induced 
measure, labeled by the size K of the environment. The probability distribution for fidelity between a random pure 
state and a random mixed state distributed according to the Bures measure is computed in Appendix B, while the 
derivation of the generating functions for the moments of the root fidelity for two mixed states of arbitrary N generated 
independently according to an induced measure is postponed to Appendix A. 



II. FIDELITY 



Let A^AT denote the set of all mixed quantum states of size TV. 
matrices of size N, which are trace normalized. 



It contains all Hermitian, semipositive definite 



Trp = 1 . 

We are going to consider two distances in this set: the Hilbert- Schmidt distance 



Dhs{pi,P2) = y/Tr[{pi - P2Y 



and the Bures distance 0, Q , 



Db{pi,P2)) = 2 - 2 Tr J ^P2^ = 2 - 2^F{pi, p^) 



(2.1) 
(2.2) 

(2.3) 



respectively. 

The Bures distance is distinguished by several remarkable properties: for pure states it agrees with the natural, 
Fubini-Study distance, while in the space of dia gon al density matrices it induces the statistical (Fisher) distance [l3l| . 
The Bures distance is a function of fidelity 0: |2^ . 



FiPl,P2) 



Tr 



P1P2VP1 



(2.4) 



This quantity is sometimes called Uhlmann transition probability Q , since for a pair of pure states it reduces to the 
squared overlap, F = Trpip2 = |(V'i|V'2)P 

Fidelity enjoys several important properties it is a symmetric, non-negative, continuous, concave function of 
both states, unitarily invariant, equal to unity if and only if both states do coincide. Therefore it becomes an important 
tool to characterize the closeness of any two mixed states, often used in modern applications of quantum mechanics 
(see e.g. ^3)- The only disadvantage consists in computing it explicitly: to find the fidelity one needs to diagonalize 
the density matrices, but more importantly, the fidelity stays being a function of the eigenvectors. 

Apart of the fidelity F, we will also use the square root fidelity y/F. Also this quantity satisfies several appealing 
properties [3 | and in certain cases its probability distribution is easier to compute then P{F). Furthermore the 
quantity arccos[-\/ F{a, p)] has a clear geometric interpretation as Bures angle in the space of mixed states equipped 
with the Bures metric. It is equal to the geodesic Riemannian distance between both mixed states 15, 16], and in 
the case of pure states it coincides with the Fubini-Study distance. 

Computation of fidelity is much simpler for N = 2. To derive an explicit formula for fidelity we use first steps of 
the calculation of the Bures distance presented by Hiibner two-by-two matrix A obeys the characteristic 

equation 



Taking the trace we obtain 



^2 _ ATtA + 1 det = 



(Tr^)2 =TtA'^ + 2dctA 



(2.5) 



(2.6) 



Let us now set A = 
of a state p. Hence 1 = 



P1P2 
(PH 



pi, so that TrA = Trpi/92 and detA 
P2f = Trp2 + 2 det p and 2 det p = 1 - 



- -y/dct pi det p2. Let (pi,p2) denote the spectrum 
Trp'^ . Substituting this into Eq. (|2.6|l we obtain 



F(pi, p2) = (TrA)2 = Trpip2 + J (l - Trp2)(l - Trp^) 



(2.7) 
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III. MEASURES IN THE SPACE OF MIXED QUANTUM STATES 



Random quantum states may be generated according to different measures, and it is hardly possible to distinguish 
the unique probability measure in the set of density matrices. Usually one considers product measures, which may 
be factorized [HIIll 

d^:E = dl^:r(Al, A2, Aat) X d/i . (3.1) 

The second factor d/i, determining the distribution of the eigenvectors of the density matrix, is the unique, unitarily 
invariant Haar measure on U {N) , while the first factor depends only on the eigenvalues Ai of p and may be chosen 
by an arbitrary probability distribution dv^ = Px{X)d^ X. The joint probability distribution is defined on the {N — 1) 
dimensional simplex Aat-i, which contains all vectors of non-negative entries summing to unity. 

The Hilbert-Schmidt measure, induced by the Hilbert-Schmidt metric, is defined by the following joint probability 
distribution [H 113 

The Bures metric leads to the Bures measure in the simplex of eigenvalues 



^N/2T{l)...T{N+l) VXiX,---Xj, y A, + A, ■ ^^-'^ 

This probability distribution was derived by Hall [Tnl |. while the normalization constants where found in |l8l | for 
N — 3, 4, 5 and in for an arbitrary N. 

Is is also instructive to consider a family of measures in the space of mixed states of size N induced by the Haar 
measure on the unitary group U{NK). The integer parameter K, used to label the measure fJ,N,K, represents the size 
of an ancilla. A mixed state of size N may be obtained by tracing a certain random pure state | $) of size NK over the 
if-dimensional ancilla, p — Tx(|$) ('I'l) • Representing the pure state |$) in an arbitrary product basis = \j) 
where i — 1, N and j — 1, K we obtain |$) = The rectangular matrix of coefficients $ij of size N x K 

allows us to write the reduced state as p — <&$'''. If matrices $ are random, the matrices p constructed in such a way 
are called Wishart matrices. The probability distribution of eigenvalues of p reads 2^ 21e 2Z | 

and ioT K = N reduces to the Hilbert-Schmidt measure H3.2() . In the above equation the inequality K > N is assumed 
and this case is called the Wishart case. If this is not the case, (the so-called anti- Wishart case 001) induced 
measures with K < N are singular, since they are supported by the subspace of states of submaximal rank K which 
belongs to the boundary of A4n- In particular, the measure fJ.N,i is just the natural Fubini-Study measure on the 
space of pure states in TYat, since the partial trace over an 1-dimensional ancilla does not change the pure state. 
Interestingly, in the one-qubit case, for K — 3/2, the induced measure H3.4|l reduces to the Bures measure (|3.3|) . 

4'^(A)-P^.,3/.(A) (3.5) 

Although such a coincidence occurs for an non-integer value of the dimensionality of the environment, and seems not 
to have any physical meaning, it will help us to obtain some results for the Bures measure by analytical continuation 
in K. Unfortunately this trick works only in the N = 2 case, for which Ai + A2 = 1, so the denominator under the 
product in the distribution (|3.3|) is trivially equal to unity. 

We are going to study mean fidelity (F) with respect to different measures. Let (F)^ denote the homogeneous case, 
in which both arguments of fidelity are random mixed states, distributed independently with respect to the same 
probability measure p. In a more general case we may average over two different measures and such averages will be 
denoted by (-F)^J. Since fidehty is symmetric F{pi,p2) = F{p2,pi), we have ~ 

An analogous notation will be also used to label probability distributions for fidelity: {F) denotes the distribution 
for the symmetric case, in contrast to Pj^^iF), in which both states are averaged over different measures. To simplify 
notation, instead of writing pn,k as a label denoting a certain induced measure, we shall use the labels K . 
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IV. MEAN FIDELITY: NON-SYMMETRIC AVERAGING 
A. One state pure, one arbitrary 

Computation of fidelity between two states simplifies considerably, if one of the states is pure, pi = \tpi){tpi\, 

J^(|V'i)(^i|,P2) -Tr|^i)(V'i|p2 . (4.1) 

Working in the basis which contains \^i) we see that fidelity is equal to the matrix element pn := {'ipi\p2\ijji) ■ We are 
going to analyze the case in which the random pure state {ipi) is generated according to the natural, Fubini-Study 
measure on the space of pure states. For any mixed state represented in a random basis all components have, on 
average, the same magnitude, hence 

(i^)FS = {F)na = ^, (4.2) 

If the second random state is also pure, p2 = \(j)){(j>\, then the fidelity is equal to the overlap between both states, 
F = \{'ip\(j))\'^ . Assuming both states to be generated independently according to the same FS measure we infer that 
fidelity is equal to the squared component of a random vector in a random basis. The probability distribution for this 
quantity 

PNAF) = {N-l){l-Ff-^, (4.3) 

was derived in [23 while analyzing eigenvectors of a random unitary matrix pertaining to circular unitary ensemble 
(CUE) and distributed according to the Haar measure on U{N). Note that this result holds also in the case if one of 
the pure states (say the second one) is fixed, P{p2) — 5{p2 — IV'*) (V'* I)- 

Let us now consider a non-symmetric averaging: one state pi = |'0i)('0i| is pure and is distributed according to the 
FS measure, while the other one, /92, is distributed according to the measure pn.k m the space of mixed states. Hence 
the latter state may be obtained by tracing a certain random pure state |$) of size KN over the X-dimensional 
ancilla. Fidelity between them, F(|-0)(?/'|,p2), is given by the sum of K terms, F — pu = J2f=i kiP where Ci, 
1 = 1, ...,KN denote the components of the pure state |$). In our previous work we have analyzed M-dimensional 
truncations of complex random vectors of dimensionality L and unit length. The length t of the truncated vector was 
shown to be distributed according to P{t) ~ ^zM-i^j^ _^2^l-m-i ^^le case considered the fidelity is just equal 

to the squared length, F = t^, the initial length of the vector L = KN, and the length of the truncation M is equal 
to K. Hence changing variables and finding the normalization constants in terms of the Euler Gamma function we 
arrive at the probability distribution 



P^-Uf) = , ^^^''"'(1 - F)^(^-i)-i (4.4) 

N.Ky ) Y{K)V[K{N ^ ' ^ ' 



describing the fidelity F between a random pure state of size TV and a random mixed state generated according to the 
measure iin,k- In the special case K — 1 the second state is also pure, and the above formula reduces to (|4.3|l . In the 
case K — N the mixed state is distributed according to the Hilbert-Schmidt measure, and the probability distribution 
reads 

For completeness, let us formulate an analogous result following from studying truncations of real random vectors 



Pnk\F) - __IiMZ!L^F^/2-l(l _ J^)K(Ar-l)/2^1_ (4 g) 

^■■^ ^ ' T{K/2)T[K{N -l)/2\ ^ ' ^ ' 



This distribution characterizes fidelity F between a reai random pure state of size N and a reai random mixed state 
generated according to the induced measure. 

Let us examine in some detail the special case of = 2. If F denotes the fidelity between p2 and than 
the fidelity of the same mixed state p2 with respect to any pure state orthogonal to is equal to 1 — F. Since we 
average over the entire Bloch sphere of pure states, the probability distribution will be a symmetric function of F 
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and (1 - F). This is indeed the case and for iV = 2 the distribution H4.4|l reduces to P^'l^ ^ ^(1 - F)^ ^ In 
particular, for the HS measure, {K = 2), and the Bures measure {K = 3/2), one obtains 
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P2,ffs - 6F(1 - F) and P2,s = (4.7) 

TT 

respectively. Eq. (|4.6|l implies analogous results for the rebits, 

PIhs = ^ and Pf,B = l. (4.8) 

Note, that for rebits P^hs and Pfg coincide, while K = 3/2 has no direct physical meaning in this case. 

In the general case of an arbitrary TV it is also possible to analyze the fidelity between a random pure state 

and a random mixed state, distributed according to the Bures measure. The probability distribution, in a sense 
complementary to Eq. 14. 5|) . 

r(jvV2)r(2iv - 1) jx - Fr'/^-^-^i - xr-y^dx 



T{N)[T{N -l/2)YT[N^/2- N) Jp ^Ny2-N+i/2 

is derived in Appendix B (to understand the method, it is more helpful first to read section V.D and appendix A). 
In the case N = 2 the integral diverges. However this divergence is compensated by the last factor, r(0), in the 
denominator, so it is not too difficult (by partial integration of the first factor under the integral) to show that H4.9|l 
reduces in this case to the second formula of (|4.7|l . A comparison of distributions of fidelity between a pure state 
and a mixed state generated according to HS measures and Bures measures is presented in Fig. ^ Observe that 
for the distributions for the Bures measure are broader since the HS measure is more concentrated in the vicinity of 
maximally mixed state. 



B. One state maximally mixed, one arbitrary 

Let us now analyze another special case, if one state is maximally mixed, pi = p, 1/N. Hence the fidelity with 
respect to any state P2 = P reduces to 

F{p^,p)^^{TT^f (4.10) 
It is then convenient to study the mean root fidelity 

=-^Trpi/2^ (4.11) 



which may be written as a function of the generalized Renyi entropy of order one half. 

Let us now assume that the random state p is distributed according to the Hilbert-Schmidt (13. 211 or Bures H3.3|l 
measure. Average moments for these ensembles of random density matrices were analyzed in [IL .12] . in which we 
derived asymptotic formulae 

(Trp^)Hs = TV^-^ . A l + Oi-)), (4.12) 



for the HS measure and 

V\CXn ^ 1 W9l / 1 \ 

(4.13) 



^ ' ' r[(i + (?)/2]r(2 + (z) V ^N' 

for the Bures measure. Substituting g = 1/2 we find asymptotic results for the fidelity of a random mixed state with 
respect to the maximally mixed state 



for the HS measure and 



^(P-P))jjg = ^ -^0.843 (4.14) 



\/2r(5/4) 



^(^-^);B^ r(3/4)r(5/2) ""-^^^ (^-^^^ 

for the Bures measure. Observe that both results, valid in the limit N — > oo, converge asymptotically to some 
nontrivial constants, independent of TV. Comparison of the Bures and the Hilbert-Schmidt measures reveals that the 
former favors states of a larger purity [TH , hence the mean value with respect to Bures measure ()4.15f) is smaller then 
the analogous result (|4.14|l for the HS measure. 
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FIG. 1: Probability distribution of fidelity between a random pure state and a random mixed state distributed according to 
the Hilbert-Schmidt measure - eq. (14.511 . solid line and Bures measure - eq. 14.91 - bold dotted line, for a) A'^ = 2, b) A'^ = 3 
and c) N = 4. 



V. MEAN FIDELITY: SYMMETRIC AVERAGING 

We are going to analyze a symmetric problem of computing average fidelity between two random states, both of 
which are generated according to the same probability distribution covering the entire set of mixed states. 

A. Case N — 2, average values 

For N — 2 the problem simplifies considerably since an explicit formula H2.7|l may be used. The expectation value 
(Trpip2) = 1/2 if pi and p2 are independent random mixed states generated with respect to any product measure 
To show this let us write both states in the Pauli matrix representation 

P^^ + r-X, (5.1) 

The vector of the normalized Pauli matrices, A = (t/\/2, together with the rescaled identity matrix, l/\/2, form 
an orthonormal basis in the set of two times two complex matrices in sense of the Hilbert-Schmidt scalar product, 
{A\B) = TiA^ B. Expanding in this basis the scalar product 



(5.2) 
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we see that due to symmetry of the distribution of the orientation of the Bloch vector f the second term does not 
contribute to the average (Tr p\P2)- For any N = 2 mixed state one has Trp^ = 1/2 — where r = \f\ denotes the 
length of the Bloch vector and r e [0, 1/V2]- Thus the fidelity may be expressed by 0| (see ea. ^2.7\i ) 



F{pi,P2) = n • f2 + ^ + - r2 - r2, (5.3) 



so to compute the mean fidelity averaged over both random mixed states distributed according to a given measure in 
the space of mixed states it is sufficient to find the expected value of the following function of the radius t — . /l — t.2 



averaged with respect to an appropriate measure. The Hilbert-Schmidt measure covers the entire Bloch ball uniformly, 
so Pffs(r) — 6\/2r^ jlflj and it is straightforward to find the average {t)HS which leads to 

{F)hs = 1 + {t)]js = i + ^vr^ « 0.6735 (5.4) 

In the same way we may evaluate the mean fidelity (F)^^ ^ with respect to the induced measures obtained by the 
partial trace of pure states of the 2K dimensional Hilbert space with an arbitrary natural K. The probability 
distribution for eigenvalues, given explicitly in implies the following radial distribution. 

Expectation values of t may be readily expressed by the Euler Gamma function and allow us to compute the average 
fidelity. For instance, in the case of the K = 3 induced measure we obtain 

{F)2,3 = l:+(^V27r) « 0.7711, (5.6) 



2 V128 

while the general result for any K > 2 reads 

1 1 / r(x + i)r(j^^i) y 

and for K = 2 reduces to (|5.4|) . In the limit of an infinitely large ancilla, X — > oo, the double quotient of the F 
functions tends to unity, so limK^oo{F)2,K = 1- This result is rather intuitive, since in this limit all random states 
tend to be close to the maximally mixed state p* . 

In the case of the Bures measure, the states of larger purity (larger radius r) are preferred and Pb (?') = ^ .[ 

[l(lllll |. Note, that the expression H5.5I) includes the above case of the Bures metric by analytic continuation: K 3/2. 
Computing the mean value (i) b we arrive at the result 

{F)b + {t)l = ^ + i « 0-590, (5.8) 

which was independently obtained by Bagan et al. and is a particular case of (|5.7|l with K = 3/2. 

It is worth to emphasize that for any product measure in the space of qubits the average fidelity is not smaller than 
1/2. The equality occurs for the Fubini-Study measure, {F)ps — 1/2, since this measure is concentrated exclusively 
on the pure states, Ppsi^) = — R), so that the last term in (|5.3|) vanishes. Here R stands for the Hilbert-Schmidt 
radius of the Bloch Ball, R = ^2/2. 



B. Case N — 2, probability distribution 



In order to calculate the probability distribution P2.k{F) we are going to use Ea. (|5.3|l and integrate out the angle 
■& between ri and T2 and the radia S [0, R] 



R rH 



P2.k{F) = ^J dn I dr2 I smM'dPK{ri)PK{r2)S 



F — rir2 cos?? — 



1 



1 



1 



(5.9) 



2 V 2 ' V 2 

The radial distribution PKif) for the measure induced by tracing over i^-dimensional environment is given in H5.5|) . 
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One first integrates out tiie angle Using the property of tiie step function 6{x — F) — d{y — F) = 9{{x — F)(F — 
y)) for X > y and performing a partial integration one finds: 

P2,k(F) = C(X)[F(1-F)]2(^-i) [ ^ (a; + 1 + 2(1 -2F))2(i-^) (5.10) 

Jq X X 

with 

CiK)^2iK-l) ■ (5.11) 

We see that the expression is analytic in K. For K ^ 1 wc obtain P2,i{F) = 1 = const. The last integral can be done 
for an integer 2{K — 1) and we obtain 



V 4 ; {2K-3)\\dFj ^F{l-F) 

Interestingly the most simple case is obtained for the Bures metric 

P2.BiF) = P2.3MF) = ^^/Fil-F) arccos(l - 2F) , (5.13) 

while for the Hilbert-Schmidt metric {K ~ 2) we obtain: 
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P2AF) = 2^(1 -F)- -^VFil-F) (1 - 2F) arccos(l - 2F) . (5.14) 

These probability distributions are presented in Fig. [3 together with exemplary numerical results. The behavior 
for K —f 00 may be obtained by a saddle-point calculation, which we start from the form 

P2,k{F) = C{K)[F{1 - ^^)]2(A'-i) i / ^ [a; + A + 2(1 - 2F)]2(i-^) (5.15) 

2 Jq X X 

using the symmetry of the integral with respect to the inversion x — ^ 1/x. We find as asymptotic form for K ^ 00 

As expected, for large dimensionality K this distribution tends to a delta function located at (F) , which in the limit 
K ^ CO tends to unity according to (|5.7|) . 



C. General case, > 3, average values 

In the general case one has to use the general formula for fidelity, (|2.4|) , not as convenient for analytical computations 
as Eq. H2.7|) valid for N = 2. Before describing the analytical results let us compare the problem for an arbitrary N 
with the symmetric case, in which the averaging in both arguments is performed over the set of pure states. For any 
nonsingular measure covering the entire set of mixed states, the mean value will be larger than 1/A'^, since the mean 
distance between random mixed states is smaller than the mean distance between random pure states. 

This statement concerns in particular all induced measures ^n,k analyzed for a fixed system size N . Since the mean 
purity decreases with the size K of the ancilla, {Tip'^) = {N + K)/{NK+1) one may infer that {F)m,k > (F) n,l 
ii K > L. In the analogy to the simplest case N — 2 one may thus expect that for any N the mean fidelity tends to 
unity, limK^oo{F)N,K = 1- 

Veracity of this reasoning may be checked by analysis of the following general results, derived in the Appendix A 
for an arbitrary induced measure ^.n,k- We are going to study the moments of the root fidelity {'\/F)n,k-, so the 
second moment gives the average fidelity. Let us denote n ^ K — N > and introduce constants 



(5.17) 



9 




FIG. 2: Probability distribution of fidelity for random one-qubit states generated according to the induced measures fJ^2,K 
where K = 1 (A), K = 2 (□), and K — 3 (0) K = 4 (o). Flat distribution represents averaging over pure states, while 
full symbols refer to the Hilbert-Schmidt measure, K = N = 2. Solid lines represent Eq. 14.31 and 15.101 . while dashed line 
denotes the case K = 3/2 corresponding to the Bures distribution 15.13(1 . 



Defining an auxiliary matrix of size N 

iXn)kj-=T{n + k + l-l)T{n + l) for k,l = 1,2,. 
the rather complicated averages may be written down in a concise way, 

{Vf)n,k = G(l) Tr[X-iX„+i/2] 

and 



N,K 



G(2) 



n+l 



Tr( (X-^X. 



n ^n+1/2) 



(5.18) 
(5.19) 
(5.20) 



The above formula is one of the main results of this paper. For N — 2 one needs to work with matrices X„ of size 
2. Computing the necessary traces one shows that (|5.2U|I simplifies to formula (|5.7|) . while (|5.19|) allows to write down 
an explicit formula 



T{2K) 



T{K) 



2\T{K -1) 



(5.21) 



For large K this average tends to unity. This result was derived for K > N, and it is ill defined, e.g ioi K — 1. 
Interestingly, by an analytical continuation — s- 1 one obtains the correct result for averaging over the space of pure 
states, {\/~F)2,i = 2/3 in agreement with the trivial integral \/FdF following from the uniform distribution l|4.3l) 
with N ^ 2. For the Bures measure we have {VF)2,3/2 = 2^13/(32527r2) « 0.7493, while the average for the HS 



measure {K = 2) is larger, {VF)2,2 = 2^^31/(5 ^7^) ft ^ 0.8098^ 

Analogous results obtained for = 3 from H5.19|l and H5.20|l read 

(VF)3,A- 



/ r(3i^) Y 


r3/r(if- 


-)\2 


3/r(if-i)x 


Vr(3if + i)/ 


L8 vr(is:- 


2)) 


A\T{K-l)) 



and 



1 



\(nK-\w 


r5 /r(if+i)x2 


1 (nK~ 


-)\2l 


Wt{k-i)) 


L12V T{K) ) 


12 vr(A:- 





2l 

2 ' • 

8 V T{K) 



i/ r(i^ + i)r(x-| )^2 

6V r{K)T{K-2) 



(5.22) 

(5.23) 
2, 



In the limit K ^ 00 both results tend to unity, e.g. {F)r^ x ^^ + "^ + 11 + ^ — 1- As in the case N 
the above formulae hold for K > N, but by analytical continuation for K ^ 1 one obtains correct results for 
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FIG. 3: Mean fidelity averaged over both states distributed according to the induced measure fiN,K, as a function of the size 
K of the ancilla ior N ^ 2 (A), N = 3 (a), N ^ A {<)), N ^ 5 (o) and = 6 (x). Sohd hues denote analytical results, Eq. 
15.711 and 15.201 1. Full symbols represent averaging over the Hilbert-Schmidt measure, K — N. 



averaging over the space of for pure states {Vf)3.i — 8/15 and (-F)3,i = 1/3. Also results in the limit K 2, 
{F)3,2 — 1/3 + 15(7r/32)^ « 0.4779 agree with the data obtained numerically by averaging over the manifold of the 
mixed states with rank K — 2. Numerical results obtained from a sample of 10^ random states, presented in Fig. 
confer with predictions obtained by means of the general formula H5.20|l . Interestingly, the mean fidelity for the HS 
measure, K = N , decreases with the system size N . 



D. General case, A'^ > 3, probability distributions 

Figure 3 presents the distributions of fidelity Pn^k{F), averaged symmetrically over both states distributed accord- 
ing to the induced measures ^n,k for iV = 2,3 and 4. Analytical results H4.3|l available for pure states {K = 1) are 
marked by solid lines. In general, the larger size K of the ancilla, the more the distribution is shifted toward higher 
values of F, since both states get closer to the maximally mixed state. For large size of the ancilla the distributions 
become concentrated at the mean value {F)]y_K which tends to unity in the limit K oo. 

Obtaining explicitly analytical results for the probability distribution P]^ x{F) in the general case of arbitrary 
A'^ and K seems not to be simple. However, we may obtain required information concerning these distributions by 
studying the sequence of higher momenta of the root fidelity. Knowing all moments {{V^)"^)n.k for to = 2,4, ... we 
may in principle extract the desired distribution Pn,k{F). To compute such moments for K > N we construct a 
generating function 



OO / A \ fJ l 

Z{\) = det(^ ^^r(^ +k + l + K-N~ l)r(^ + 1 + K-N)) (5.24) 

m— 



which we obtain by replacing the fixed trace ensembles for p\ and pi by Laguerre ensembles and integrating out 
eigenvectors of ^/p\Pi^fp\ with the help of the Itzykson-Zuber integral. The determinant concerns a N x N matrix 
with indices, k,l = 1, .., N. To obtain the mean values of (VF)™ one needs to multiply the expansion coefficients of 
Z{X) by constants G{m) defined in H5.17|l : 



{WFy^)N,K = G(to) i-d/dXr Z(A)/Z(0)U=o (5.25) 

For any fixed N the moments are analytic in K, so one may try to use the analytic extension for the cases K < N. 
Alternatively, we may treat this case separately, investigating another generating function 

Z.iX) = detf V [r(™ +k + N-K)r + ^ + (5.26) 

' ml L V2 )i r('4+k + l + N-2K)J ^ ' 



11 




0.5 F 1.0 




0.5 F 10 



FIG. 4; Probability distribution of fidelity while averaging symmetrically over both states distributed according to induced 
measures fiN,K- The size of the system equals a) A'^ = 3 and b) = 4 and while the size K of the ancilla is coded by symbols: 
K = 1 (A), K — 2 (□), K = 3 (0) and K = 4 (o). Full symbols represent averaging over Hilbert-Schmidt measure, K = N. 
Solid line for 7^ = 1 denote Eq. 114.31 , dotted lines are plotted to guide the eye. 



with determinant oi a. K x K matrix with indices, k,l — 1, .., K. This formula holds for N > 2K, but may be extended 
analytically in N beyond this restriction. It is obtained by reducing ~/p\p2\fp\ by unitary rotations to the spaces of 
nonzero eigenvalues with dimension K < N using the symmetry properties of the fidelity and again integrating out 
all unitary rotations with the help of the Itzykson-Zuber integral, (here applied two times). 

To show an application of this approach let us expand the generating function up to the second order in A. 
Fixing K = 2 we obtain the averages for variable iV 



^)Na = ^ (22iV-13) ( 



V{2N) N^2 r(iv- i) 

r(iv) 



r 27V- 



and 



1 r 



{F)n.2 

These formulae are valid for arbitrary iV, and for N 



1 



Stt T(N 
16 



\)T{N- 



(5.27) 



(5.28) 



r(iv- i)r(7v + 1) 

2, 3 are consistent with Eqs. H5.21|l - H5.23|l . In a similar way one 
may proceed with K = 3,4.... Furthermore, expanding the generating function to the m-th order one could calculate 
the expectation values of (-s/F)™ and obtain further information concerning the probability distribution Pff^KiF). 
There are analogous formulas for the moments ((a/F)™) n,k to Eqs. (|5.19|l . 15.20|l in terms of traces involving K x K 
matrices. 

To give an explicit expression for the distribution of fidelity Pn^k{F) = W{VF)/2VF we start from the fixed trace 
ensemble with traces ti = Trpi, t2 — Trp2, ^3 = '^^{y/PiP2\/piy^'^ ■ We take the Laplace transform with respect to 
ill ^2 7^3 J calculate it and then apply the inverse Laplace transform at ti = 1, t2 = 1, and = \/F. What we obtain 
is a threefold contour integral 



W{VF) ^ C 



dsie'*! 
27risi^^ 



ds2e'^ 
2ms2^^ 



2TTi 



S3 



y/sTs2 



(5.29) 



with C = r{K'^N'^)/Z{0). It is easy to see, that this equation is equivalent to Eq. 15.25|l . All three contours go along 
the whole imaginary axis with a small positive real part. We need the analytic properties of the function Z(x) which 
can be seen from the form 



Zix) 



y 



fe-2 



T[k-l + 2iK - N + 1)] 



det ( 2 1^ dy- ^ ^ ^Y~^+2(K-N+i) 



(5.30) 



Z{x) has a cut along the real axis for x < —2 and behaves asymptotically for x — > cx) like (lnx)^/x^^^, which is not 
easy to see, but can be derived from an explicit form for Z{x) 



Z{x) = det[(d/dx)2(^-^+'-i)+^-i/fc(x) ln(a;/2 + ^2^2/4-1)] 



(5.31) 
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with 

Hk{x) = ~ [{x/2 + ^xy^-l)'"'' + {x/2 - v/^V4^)fe-i] /^a;2/4_i (5.32) 
For later use we need also for a; > 2 (in the following e > 0) 

lim[Z{-x-ie)] =det[(d/da;)2(^-^+'-i)+'=i/fc(2;)ln(a;/2+ v/a;2/4- 1) + nT{d/dxf^^-^+^-^^+^ Hk{x)] . (5.33) 

Using appropriate contour deformations and partial integrations we are able to perform the si and S2 integrations 
arriving at: 



with 



2/%/F 

KN~l ' 



B(x) = F 



dxB{x)Im[Z(—x — ie)] 
dy e{x-y) {x-yf^^-^ 



(5.34) 



(5.35) 



Finally for numerical convenience we may use the following integral representations for x > 2: Z{x) — det(Ai) and 
lini£^o[^(~2; — ie)] — det(Ai — 1^2) with 



and 



'° du cosh[(/c - l)u]T[k -l + 2iK -N + 1)] 
[2cosh(ii) +a;]'=-i+2(^<'-^+0 

" ducos[{k - l)u]r[k - 1 + 2{K -N + 



[2 cos(w) + a;]fe-i+2(-^f-^^+0 
where k,l = \, . . . , so A\ and A2 are square matrices of size N . 



(5.36) 



(5.37) 



VI. CONCLUDING REMARKS 



In this work we have posed and solved the problem of computing the average fidelity between two random quantum 
states. As easy to predict, the results depend heavily on the probability measure, according to which the quantum 
states are generated. We have concentrated on two cases, we consider to be the most important: the Bures measure 
/iB, used to guess a random state of size N in lack of any other information, and the induced measure fJ-N,K, applied 
if it is known that the mixed state is obtained by partial trace over a K dimensional environment. The latter case 
reduces to the Hilbert-Schmidt measure if K = N. 

In general, one has to distinguish the case of the symmetric averaging, in which both states are generated according 
to the same measure, and the asymmetric case, e.g. the problem of computing the mean fidelity between a pure 
random state and a mixed random state. 

In particular we have obtained explicit results for symmetric averaging over the induced measures fJ.N,K, which have 
a simple physical interpretation. The mean value (F) 7v,if corresponds to the case of picking at random pure states 
in iV X if dimensional Hilbert space and then studying the fidelity between both mixed states obtained by partial 
tracing over the if-dimensional environment. 

In certain cases we have derived the probability distribution for fidelity P{F) which may be used to formulate 
statistical statements which evaluate quantitatively certain quantum operations and quantum protocols. For instance, 
the universal protocol of quantum cloning of a dimensional pure state gives the fidelity Fu = {N + 3)/{2N + 2) 
|28j . Making use of the distribution (|4.3|) we see that the probability pn that a randomly taken N dimensional pure 
state will give a better fidelity than this obtained by the cloning procedure is equal 

1 /A-l\^-i 

P--Jj-AP)^P-W-^{j^) • (6-1) 

For N — 2 there exists a considerable probability, p2 = 1/6, that picking at random a pure state of a qubit provides 
results better then the cloning protocol. However, for larger system sizes the probability decreases exponentially 
with N , which means that the relative strength of the cloning procedure with respect to the 'random choice' strategy, 
increases with the size of the system. 

We thank D. Savin for useful discussions and I. Bengtsson, R. Mufioz-Tapia and A. Uhlniann for helpful correspon- 
dence. This work was supported by Sonderforschungsbereich / Transregio 12 der Deutschen Forschungsgemeinschaft 
and a solicitated grant number PBZ-MIN-008/P03/2003 of Polish Committee for Scientific Research. 
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APPENDIX A: GENERATING FUNCTIONS FOR MOMENTS 
1. Induced measures: Wishart case, K > N 

In this appendix we present a brief derivation of the generating functions for moments of the square-root fidehty 
which hold for an arbitrary induced measure ^n,k- In the first subsection we consider the case K > N. Then the 
distribution of square-root fidehty is given by 

WiVF) J Dpi J Dp2 det{pip2f-''S{l - Trpi) 6(1 ~ Trpa) s(^Vf ^ Tv{^p2^)^^^) , (Al) 

where Dpi, Dp2 are the matrix volume elements and we have to integrate over positive matrices pi and p2- Since 
Tr{{yfpiP2y/pi)^^'^) is a homogeneous function of both, pi and p2 of degree 1/2, moments can equivalently be obtained 
from the generating function 

Z{\) J Dpi J Dp2 dctipip2f'''e-^'P^ e-^P- e-^^r[i^P2VpT)'^'] (A2) 

following Equ. (|5.25|) . Wc make a transformation p2 Pi^^^^ P2Pi^^^^ and taking into account the Jacobian we 
obtain 

Z{X) oc Jdpi J Dp2 dct(pi)-^dct(p2)^~'^e-T''^i q-^'^p^/pi e-ATr(p.)^/^ (A3) 

Note that the nasty square roots disappeared and we can use the Itzykson-Zuber integral to integrate out the 
eigenvectors of p2 in the expression q-^'^p^/pi. As result we obtain 



Z{X) cx 



/ dyi...dyiv(2/i...y^)^-^e-(-^+--~)e-^(v17r+-v^) detie-y^^^^')ri{x,-xi)iy,~yi) (A4) 
Jo XI xr, J 



where xi, ...xn are the eigenvalues of pi and yi, ...yN are the eigenvalues of p2- Using the properties of the Vander- 
monde determinant we reduce this expression to a single determinant and arrive at Eq. H5.3U|) . 

2. Induced measures: anti— Wishart case, K < N 

In this case the matrices pi and p2 are of rank K and they can be brought by a unitary transformation of the form 

\ Vi-x^xJ 

to block diagonal form. If U is drawn from CUE it can be shown that the corresponding measure DX is flat, where 
X is a complex K x {N — K) matrix with the restriction XX'^ < 1. Calling the nonzero blocks again pi and p2, now 
K X K matrices, we may write the generating function for moments as 

Zk{X) oc [dpi I Dp2 I DX det(pip2)^-^e-T''^^ e-^'''^ ^-XTr[i^,^i-xxi p.^i-xx^ ^Y'-] (^5) 

In a first step we make the transformation P2 ^ (1 ^ XX^)^^/^p2{l — XX^)~^^'^ arriving at 

Zk{X) J Dpi J Dp2 j^^^^DX ^^ElP^^o-^^P^ ^-T.p./ii-xxr) ^-ATrKvT^Tp.v^)^/^] . (A6) 

Then we use the symmetry ^^^:{{^/plP2^/PlY^^) — Tr((y^piy^)^/^) and make the transformation pi 
P2^^^'^ P\P2^^^'^ leading to 
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Now we see that the matrix W = XX^ is again of the Wishart type and is positive definite for K < (TV — K). Let us 
focus on this case, then the distribution of is oc det(W^^~^^) and we may reduce the integral to the eigenvalues 
xi, ...xk of Pi, yi, ...yK of p2 and zi, ...zk of W using two times the Itzykson-Zuber integral. The resulting expression 
contains in the integrand the Vandermonde determinants of xi^...xk and of zi^...zk and can therefore again be 
reduced to a single determinant which is given by 



Zk{X) = det 2 



dx 



dy 

y 



dz- 



(1-z) 



N-K-l+l 



■ exp 



X f 1 

— {y + - 



(A8) 



Expanding the threefold integral in powers of A, one arrives at Eq. (|5 



APPENDIX B: DERIVATION OF THE DISTRIBUTION 

To derive the distribution P{F) for the fidelity between a random pure state and a random mixed state distributed 
according to the Bures measure 1)3. 3|l consider 

P{F)oc [ Dp [ DA S{l-TTpA^) S{l-Trp) S{F -TTp\iP){ij\), (Bl) 
Jp>0 Ja=a^ 

where we introduced an auxiliary Hermitian matrix A to generate the distribution 1)3. 3|1 and the last two (5-functions 
correspond to the constraints H2.1|l and (|4.1|) . The J- functions may be written as inverse Laplace transforms along 
the slightly to the right shifted imaginary axis and then the integration over the positive definite iV x A^-matrix p can 
be done arriving at 

After a suitable rescaling of the matrix A, we can do the si- and sa-integrations and obtain 

» / d.t(iT.4.)~ a ^^^^^ ^ ■^''"'"»'" ""'^ ''•K-whSfm) ^ 

Now it is possible to show, that with the correspondence U — {1 + iA)/{l — iA) the measure _D^/det(l + A^)^ 
is equivalent to the invariant measure on the group of unitary A^ x A-matrices U. Since only {tp\{l + A'^)'^\4>) = 
{ip\{2 + U + U^)/4:\tp) occurs, the integration over A can be reduced to a one-dimensional one and after integrating 
over s in the complex plane and restoring the normalization we end up with formula H4.9|l . 
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